UAB-FT-653 



Ktt vector form factor, dispersive constraints 

§ and T -^ Ut-Ktz decays 

o' 

(N 

o 

^ \ Diogo R. Boito^, Rafel Escribano^, and Matthias Jamin^'^ 

(N 



o 



^ Grup de Fisica Tedrica and IFAE, Universitat Autonoma de Barcelona, 
E-08193 BcUatcrra (Barcelona) , Spain. 



Oh' ^ Institucio Catalana dc Rcccrca i Estudis Avangats (ICREA). 

(D ' 

£^ 

(N 
> 

QQ ' Abstract: Recent experimental data for the differential decay distribution of the de- 

OO ■ cay T~ — > I'rKsTT' by the Belle collaboration are described by a theoretical model which 

is composed of the contributing vector and scalar form factors F:^'"{s) and Fq^'^{s). Both 

form factors are constructed such that they fulfil constraints posed by analyticity and unitar- 



OO , ity. A good description of the experimental measurement is achieved by incorporating two 

vector resonances and working with a three-times subtracted dispersion relation in order to 
suppress higher-energy contributions. The resonance parameters of the charged K*{892) me- 

r> I son, defined as the pole of F^'^(s) in the complex s-plane, can be extracted, with the result 

C^ ■ Mk* = 892.0 ± 0.9 MeV and Tk* = 46.2 ± 0.4 MeV. Finally, employing a three-subtracted dis- 

persion relation allows to determine the slope and curvature parameters A^ = (24.7ib0.8)-10~'^ 
and A^ = (12.0 it 0.2) • 10~^ of the vector form factor F:^'^{s) directly from the data. 



PACS: 13.35.Dx, ll.30.Rd, 11.55.Fv 

Keywords: Decays of taus, chiral symmetries, dispersion relations 



1 Introduction 

Hadronic decays of the r lepton provide a fruitful environment to study low-energy QCD 
under rather clean conditions [1-5]. Many fundamental QCD parameters can be determined 
from investigations of the r hadronic width as well as invariant mass distributions. A prime 
example in this respect is the QCD coupling a^ [6-8]. In addition, fundamental parame- 
ters of the strange sector in the Standard Model can also be obtained from the r strange 
spectral function. The experimental separation of the Cabibbo-allowed decays and Cabibbo- 
suppressed modes into strange particles [9-11] paved the way for a new determination of the 
quark-mixing matrix element |V^s| [12-15] as well as the mass of the strange quark [16-23]. 

Cabibbo-suppressed r-decays into strange final states are dominated by r — > Vt-Ktt. In 
the past, ALEPH [9] and OPAL [10] have measured the corresponding distribution function 
but lately i?-factories have become a new source of high-statistics data for this reaction. 
Recently, the Belle experiment published results for the r -^ Vt-Ktt spectrum [24] and a new 
determination of the total branching fraction became available from BaBar [25-27]. In the 
future, there are good prospects for results on the full spectrum both from BaBar and BESIII. 

Theoretically, the general expression for the differential decay distribution of the decay 
T — > Vt-Ktt can be written as [28] 



V^ - 32vr3s ^^^V M^j 



l + 2^\\^\F^-{s)\^ + ^-^qK.\F,^-{s)? 



(1) 
where isospin invariance is assumed and we have summed over the two possible decay channels 

T~ — > VrK TT~ and r~ — > VrK~T^^ , with the individual decays contributing in the ratio 2 : 1 

respectively. Furthermore, Sew = 1-0201 [29] represents an electro- weak correction factor, 

Ai^TT = M^ — M^ , and qxir is the kaon momentum in the rest frame of the hadronic system, 



QkAs) = ^Y (^ - i^^K + M^f) [s - (Mk - M^)2) . e(^s - {Mk + M^f) . (2) 

Finally, we denote by F^'^{s) and Fq'^{s) the vector and scalar K-k form factors respectively, 
which we will discuss in detail below. 

In eq. ([1]), the prevailing contribution is due to the K-k vector form factor Fi^'"{s), and in 
the energy region of interest, this form factor is by far dominated by the i^*(892) meson. A 
description of F^'^{s) based on the chiral theory with resonances (RxT) [30,31] was provided 
in ref. [32], analogous to a similar description of the pion form factor presented in refs. [33-35]. 
Then in ref. [36] this description was employed in fitting Belle data for the spectrum of the 
decay r — > z/t-X^tt" [24]. The additionally required scalar Ktt form factor Fq'^{s) had been 
calculated in the same RxT plus dispersive constraint framework in a series of articles [37-39] , 
and the recent update of F^'^{s) [40] was incorporated as well. 

A slight drawback of the description for the vector form factors of refs. [32, 33] is that 
the form factors only satisfy the analyticity constraints in a perturbative sense, that is up 
to higher orders in the chiral expansion. Though the violation of analyticity is expected to 
only be a small correction (of order p^ in the chiral expansion in the case at hand) it is 



certainly worthwhile to corroborate this assumption. A coupled channel analysis of the Kir 
vector form factor, which would allow for such a test, was performed in ref. [41]. However, 
in ref. [41] the theoretical description was not really fitted to the experimental data, so that 
it is difficult to decide if differences of the coupled channel analysis [41] as compared to the 
description [32, 36] already show up in the current experimental data. Even more so as the 
fits performed in ref. [36] provided a satisfactory description of the Belle spectrum. 

Below, we shall investigate the related questions in a more modest approach. In the region 
of the K*(892) meson, elastic unitarity is still expected to hold. Since this meson dominates 
the Kir vector form factor, an ansatz implementing elastic unitarity should result in a good 
approximation. For the pion vector form factor such an approach was pursued in ref. [34] 
and in the present work we perform an analogous investigation for F:^'^{s). Even though 
possible coupled-channel contributions are not explicitly included in our parametrisation of 
the Kit vector form factor, their influence can be studied through the sensitivity of our ansatz 
when changing the number of subtractions in the dispersion relation that the form factor 
satisfies, because a larger number of subtractions entails a stronger suppression of higher- 
energy contributions. As a benefit of our approach, we are able to extract the resonance 
parameters of the K*{892) from the pole of F:^'^{s) in the complex s-plane, which should be 
regarded as more model-independent than the Breit-Wigner type parameters extracted in the 
previous analyses [24,36], as well as the first three slopes in the Taylor expansion of F^'^{s) 
around s = 0. 

2 The Ktz vector form factor 

The Kit vector form factor F:^'^{s) is an analytic function in the complex s-plane, except 
for a cut along the positive real axis, starting at the Kn threshold skw = {Mk + Mt^)'^, where 
its imaginary part develops a discontinuity. The analyticity and unitarity properties of the 
form factor result in the fact that it satisfies an n-subtracted dispersion relation, explicated 
in more detail for example in refs. [33,34]. In the elastic region below roughly 1.2 GeV, the 
dispersion relation admits the well-known Omnes solution [42] 

^ oo -. 

F^^s) = Pn{s) expi ^ fds' ,^,J',^^''} ,,, , (3) 



vr J (s')"(s'-s-iO) 



which corresponds to performing the n subtractions at s = 0, and where 

I (4) 



I fc=0 



is the subtraction polynomial. More general formulae with subtractions at an arbitrary point 
s = sq can for example be found in ref. [43]. Furthermore, Si{s) is the P-wave / = 1/2 
elastic Kir phase shift. As on general grounds the phase shift is expected to go to an integer 
multiple of vr, at least one subtraction (n = 1) is required in eq. ^ to make the integral 
convergent. Let us first dwell on this case in more detail, before turning to the case with a 
larger number of subtractions. 



While in the approach of refs. [32, 33] to the vector form factor the real part of the one- 
loop integral function H{s) (for its definition see [32]) is resumnied into an exponential, strict 
unitarity is only maintained if this piece, together with the imaginary part which provides 
the width of the resonance, is resummed in the denominator of the form factor [37,44]. The 
resulting expression of the vector form factor corresponding to one single resonance then reads 

F^^s) = - <* (5) 

m^. - s - kHktt{s) 

Of course, up to order p^ in the chiral expansion, resumming the real part in the denominator 
or an exponential is fully equivalent. Differences of both approaches first start to appear at 
0{p^). In eq. ([5]) the parameter rrix* is to be distinguished from the true mass of the K* 
meson M^*, which later will be identified with the real part of the pole position of F^'^{s) 
in the complex s-plane|l| 

Identifying the imaginary part in the denominator of eq. ^ with —mK*jK*is), where 
the s-dependent width of the K* meson takes the generic form of a vector resonance, 

with 7i^. = 7i<-.(m,|^*), the dimensionful constant k has to take the value: 

1927rFKF^ jK* 



(TKniml^*)^ rriK* 



(7) 



In eqs. ([6]) and ([7]), the phase space function crKiris) is given by axnis) = 2qKTT{s)/ \/s. The 
form factor F^'^{s) can therefore be written in the equivalent form 

F^^{s) = ^ . (8) 

m?i^t — s — KKeHKnis) — i'mK*lK*{s) 

From eq. dH]) the normalisation F_p'^(0) at s = which is needed in order to calculate the 
reduced form factor Fi^'^{s) = F^'^{s)/F:^'^{0), is given by 

F^^O) = ^^^^h-— (9) 



The phase 5i{s) of the form factor F^'^{s) is found to be: 



,K.^ ^ _ ImFf'-(s) niK'^K^is) 



tan(5f'^(s) = ±--J- = ^ '" ^ 1 . (10) 

^ "■ ' ReFf-(s) m\,-s-KReHK^{s) ^ ' 

It is a straightforward matter to verify that with the phase 6i{s) of eq. (llOp . the form factor 
of eq. ([8]) satisfies the Omnes relation ^ with n = 1 and -Pi(s) = Ff^'^{0). Therefore, it is 



^The renormalisation scale fj, appearing in HKn{s) will be set to the physical resonance mass fi = Mk* ■ 



in accord with the analyticity and unitarity requirements. Finally, the reduced form factor is 
written as _ 

F¥'(s) = %^4f<5^ . (U) 

nij^, — s — KtieHKn{s) — imK*^K*\S) 

which resembles the non-strange form factor of Gounaris and Sakurai [45]. (For comparison, 
see also ref. [46].) However, while in the Breit-Wigner resonance shape of ref. [45] the real part 
of the inverse propagator is obtained through a twice subtracted dispersion relation combined 
with proper subtractions to fix mass and width of the resonance, the form factor in eq. (jlip 
is found from a once subtracted dispersion relation satisfying analyticity and unitarity. 

As has already been discussed in ref. [36], in eq. ([T]) only the reduced form factor F^'^{s) 
has to be modelled, as the normalisation oi Fi^'^{s) only appears in the product |Kjs|-P+'^(0). 
This combination is determined most precisely from the analysis of semi-leptonic kaon decays. 
A recent average was presented by the FLAVIAnet kaon working group, and reads [47] 

iKsli^+^^^lO) = 0.21664 ±0.00048. (12) 

In what follows, we work with form factors normalised to one at the origin and assume the 
value (|12p for the overall normalisation. We remark that the normalisation for the scalar and 
the vector form factors is the same and that (J12p already corresponds to the K^tt~ channel, 
which was analysed by the Belle collaboration [24]. Consequently, possible isospin-breaking 
corrections to F^ '^ (0) are properly taken into account. 

3 Single resonance fits to the Belle spectrum 

Our fits to the Belle r" — > u^Kst^^ spectrum [24] will be performed in complete analogy to 
the recent analysis of ref. [36]. Let us briefly review the main strategy for these fits. The 
central fit function is taken to have the form 

\-\- 0.0115 [GeV/bin] AT^r • :-i— ^ • (13) 

2 3 TrBxTT d^s 

The factors 1/2 and 2/3 arise because the Kstt~ channel has been analysed. Then, 11.5 MeV 
was the bin-width chosen by the Belle collaboration, and At = 53110 is the total number of 
observed signal events. Finally, F,- denotes the total decay width of the r lepton and Bk-k a 
remaining normalisation factor that will be deduced from the fits. The normalisation of our 
ansatz (J13p is taken such that for a perfect agreement between data and fit function, Bxn 
would correspond to the total branching fraction Bxn = B[t~ -^ z/^A's'Tt"] which is obtained 
by integrating the decay spectrum. All further numerical input parameters have been chosen 
as in ref. [36]. 

Numerically, our first fit of the Belle data [24] with the vector form factor F^'^[s) according 
to eq. dS]), and including data up to ^/s = 1.01525 GeV (centre of bin 34), is presented in 
tabled! As discussed in ref. [36], we have removed the problematic data points 5, 6 and 7 
from the fit. Furthermore, we have not taken into account the lowest data point since for our 
physical meson masses the centre of this bin lies below the Kir threshold. For the scalar form 
factor Fq'^{s) we have used the recent update [40], employing the central parameters given 



Eq. dlD for F_f ^(s) 



Bktt {Bktt) 0.3611 ± 0.0042 % (0.3562 %) 
niK* 943.35 ± 0.51 MeV 

7A-* 66.29 ±0.79 MeV 

xVn.d.f. 39.4/27 



Table 1: Results for the fit with a one-subtracted dispersion relation including a single vector 
resonance in Fj^'^{s) according to eq. ([8]), as well as the scalar form factor Fq'^{s) [40]. 



there. One observes that due to the not very satisfactory quality of the fit, the fit parameter 
Bktt and the integrated branching fraction Bxn display a marked deviation, which is however 
within the uncertainties. Furthermore, because of the real part of the loop-integral in the 
denominator, mx* and 'Jk* turn out rather different from their physical values Mk* and 
Tk*- For the final results of the parameters in our description of the vector form factor, 
in our conclusion we shall also present values for the physical parameters Mk* and Tk*, as 
obtained from the pole position in the complex s-plane. 

In order to make the fit less sensitive to deficiencies of our description in the higher-energy 
region, a larger number of subtractions can be applied to the dispersion relation. Besides, 
employing an n-subtracted dispersion relation has the advantage that the slope parameters 
which appear as subtraction constants are determined more directly from the data. It should 
be pointed out, however, that the form factor with a larger number of subtractions n > 2 
violates the expected QCD large-energy behaviour. For large enough energies the form factor 
F^'^{s) should vanish as 1/s. Since we only employ the vector form factor up to about 
^/s ~ 1.7 GeV, which is still in the resonance region, we consider this deficiency acceptable. 
Anyhow, as can be verified explicitly from our fits, in the considered region above the second 
vector resonance, F^'^{s) is a decreasing function of s. On the other hand, fits with only one 
subtraction were generally found to only provide a poor description of the experimental data. 

Let us present our results for the case n = 3 in detail, but below, we shall also briefly 
comment on the cases n = 2 and n = 4. The importance of the high-energy region can be 
studied by introducing a cutoff Scut as the upper limit of the integration in the Omnes integral 
([3]). Incorporating three subtractions, the reduced form factor F^'^{s) then takes the form: 



(14) 
I ivj "_ z ivi '_ 7( J \^ti- j-'^^- — if — lyj) I 

where the phase 5f'^{s) corresponds to the expression of eq. p^ . The parameters ai and a2 
can easily be related to the slope parameters A^_ , which appear in the Taylor expansion of 
F^'^{s) around s = : 

— T^^ , -. f S 1 // S 1 '" s , , 

F^^(s) = 1 + X,^^ + -X,^^ + -X,^^ + ... . (15) 

+ V ; -r + j^2_ ^ 2 + M^_ 6 + M^_ ^ ' 




Explicitly, the relations for the linear and quadratic slope parameters A_)_ and A_)_ then take 
the form: 

A_|_ = ai , A_|_ = 02 + «! • (16) 

Below, we shall also compute the cubic slope parameter A, from the dispersive integral. 





Scut = 3.24 GeV^ 


seut = 4 GeV2 


Scut = 9 GeV2 


Scut -^ OO 


Bk. 


0.394 ± 0.045 % 


0.397 ± 0.046 % 


0.398 ± 0.046 % 


0.398 ±0.046% 


{Bk.) 


(0.389 %) 


(0.391 %) 


(0.393 %) 


(0.393 %) 


rriK* [MeV] 


943.34 ± 0.57 


943.36 ± 0.58 


943.37 ±0.58 


943.37 ± 0.58 


7i^. [MeV] 


66.48 ± 0.87 


66.50 ± 0.89 


66.52 ±0.89 


66.52 ±0.89 


a; X 103 


23.9 ±2.4 


24.2 ±2.5 


24.4 ±2.5 


24.4 ± 2.5 


a; X 10^ 


11.5 ±0.7 


11.5 ±0.7 


11.5 ±0.7 


11.5 ±0.7 


xVn.d.f. 


45.8/41 


45.8/41 


45.7/41 


45.8/41 



Table 2: Results for the fits with a three-subtracted dispersion relation including a single 
vector resonance in Fi^'^{s) according to eq. ([8]), as well as the scalar form factor Fq'^{s) [40]. 

The results of our fits with the three-subtracted dispersion relation, employing four values 
of i/Scut , namely 1.8 GeV, 2 GeV, 3 GeV and Scut -^ oo, are given in table [2j For these fits, 
we have included experimental data up to the data point 50 at ^/s = 1.19925 GeV (centre of 
the bin) [24], and as already mentioned above, removing the problematic data points 5, 6 and 
7. As can be seen from table O the fits are indeed rather insensitive to the upper integration 
limit Scut, implying that the higher-energy region is well suppressed. Compared with the fit 
of table [H which only employed a single subtraction, with x^/n.d.f. ~ 1.1 the fit quality is 
substantially improved. Related to that, also the results for Bk-k and Bk-k turn out much 
closer. However, the slope parameter A^ is not well determined from our fit. This originates 
in the fact that the slope parameters are almost 100 % correlated with the total branching 
fraction Bk-k-, and this parameter has relatively large uncertainties. Therefore, in our best 
estimate of the model parameters below, we shall impose the experimental measurement of 
the total branching fraction Bktt- 

Though we only present explicit results for the three-subtracted dispersion relation, we 
have also investigated the cases n = 2 and n = 4. In the case n = 2, a still somewhat stronger 
dependence on the cutoff Scut is observed, which is why we do not discuss the corresponding 
results in detail. The four-subtracted dispersion relation, on the other hand, yields almost 
unchanged central values for the fit parameters without any improvement in the x^/n.d.f., but 
with larger parameter errors due to the additional degree of freedom. From this we conclude 
that the case n = 3 discussed above is an optimal choice as far as the number of subtractions 
is concerned. Next, we shall also include a second vector resonance, the i('*(1410) into our 
description of the Ktt vector form factor. 



4 Fits with two vector resonances 

For our final fits, we aim at a description of tlie r — > v-j-Kir spectrum in the full energy 
range up to the r mass. To this end, we also introduce a second vector resonance, the 
K* = if* (1410), into our model for the vector form factor F:^'^{s). As it is unclear how 
to directly incorporate a second vector resonance in the phase 6i{s), we have followed a 
somewhat indirect approach, which should be sufficient for our purposes, as the contribution 
of the K* resonance is suppressed by phase space. 

Like for the case of the single resonance, as a starting point we assume the form of F^'^[s) 
given in eq. (5) of ref. [36], however resumming the real part of Hk-k{s) in the denominator 
of the form factor. This leads to the following expression: 

where 

D{mn,jn) = m^- s - KnReHKn{s) -irrinjnis) ■ (18) 

For both resonances, "fnis) is given equivalently to the form of eq. ([6]), and the corresponding 

K„ can be deduced in analogy to eq. ([7]). Like in eq. ([TO]) , the phase 5^'^{s) can be calculated 

from the relation 

^r.^, , ImFi^''^(s) , , 

This is the phase that we then employ in the Omnes integral representation ([3]) for the form 
factor to perform our fits. 





s,„t = 3.24 GeV^ 


Scut = 4 GeV^ 


Scut = 9 GeV2 


Scut -^ OO 


Bk^ 


0.386 ± 0.043% 


0.404 ± 0.044% 


0.417 ± 0.046% 


0.417 ± 0.046% 


(Bk.) 


(0.384%) 


(0.402%) 


(0.414%) 


(0.414%) 


rriK* [MeV] 


943.27 ±0.58 


943.40 ± 0.57 


943.48 ± 0.57 


943.49 ± 0.57 


7i^. [MeV] 


66.43 ± 0.90 


66.66 ±0.87 


66.81 ±0.86 


66.81 ±0.87 


m^.' [MeV] 


1392 ± 41 


1369 ± 30 


1361 ± 28 


1361 ± 28 


iK-' [MeV] 


273 ± 137 


224 ±101 


212 ± 93 


212 ± 93 


7X 10^ 


-4.2 ±2.0 


-3.6 ±1.6 


-3.4 ±1.5 


-3.4 ±1.5 


A'+ X 103 


22.6 ±2.2 


23.9 ±2.1 


24.7 ±2.1 


24.8 ±2.1 


a; X 10^ 


11.5 ±0.6 


11.7 ±0.7 


11.9 ±0.7 


11.9 ±0.7 


xVn.d.f. 


73.7/78 


75.6/78 


77.2/78 


77.3/78 



Table 3: Results for the fits with a three-subtracted dispersion relation including two vector 
resonances in F^'^{s) according to eqs. ()17p to ()19p . as well as the scalar form factor Fq'^{s) 
of ref. [40]. 

Our first fit with two vector resonances proceeds in complete analogy to the second fit 
in the single-resonance case. Again, we use the three-subtracted dispersion relation, and 



investigate four values of Scut, in order to study the dependence of our fits on the higher- 
energy contributions. Now, the Belle data [24] have been included up to ^/s = 1.65925 GeV 
(data point 90) o The corresponding fit results are presented in table [3l 

The general picture of our fits with two resonances is quite satisfying. The theoretical 
model provides a good description of the experimental data in the full energy range. The 
X^/n.d.f. for all values of Scut turns out smaller than one. The dependence of the resulting 
fit parameters on Scut is small and within the uncertainties, though clearly visible. The fits 
provide a precise determination of the parameters of the lowest lying K* vector resonance, 
and a still reasonable accuracy for the second K* resonance. However, the fit uncertainties 
for the branching fraction Bktt are relatively large. Accordingly, also the error on the slope 
parameter A_|_ is found large, because these two parameters are almost 100 % correlated. 

Scut = 3.24 GeV2 s^ut = 4 GeV^ s^nt = 9 GeV^ Scut ^ oo 



rriK* 


[MeV] 


943.32 ± 0.59 


943.41 ± 0.58 


943.48 ± 0.57 


943.49 ± 0.57 


IK* 


[MeV] 


66.61 ± 0.88 


66.72 ± 0.86 


66.82 ± 0.85 


66.82 ± 0.85 


ruj^. 


/ [MeV] 


1407 ± 44 


1374 ± 30 


1362 ± 26 


1362 ± 26 


Ik*' 


[MeV] 


325 ± 149 


240 ± 100 


216 ±86 


215 ±86 


!'■ 


K102 


-5.2 ±2.0 


-3.9 ±1.5 


-3.5 ±1.3 


-3.5 ±1.3 


a; 


X 10^ 


24.31 ± 0.74 


24.66 ± 0.69 


24.94 ± 0.68 


24.96 ± 0.67 


Al 


X 10^ 


12.04 ±0.20 


11.99 ±0.19 


11.96 ±0.19 


11.96 ±0.19 



xVn.d.f. 74.2/79 75.7/79 77.2/79 77.3/79 

Table 4: Results for the fits with a three-subtracted dispersion relation including two vector 
resonances in F^'^{s) according to eqs. p!7j) to p^ . as well as the scalar form factor F^'^[s) 
[40]. The total branching fraction Bk-k has been fixed to the experimental value (I20p . and 
the corresponding uncertainty is included in quadrature. 



In order to provide a better determination of the slope parameters, a possible way to 
proceed is to fix the total branching fraction Bk-k to the experimental measurement. A very 
recent update of the world average has been presented in ref. [24,26,27], with the finding: 

B[t- -^ VrKsTT~] = 0.418 ± 0.011 % . (20) 

Thus, for our final fits, we have fixed Bk-k to take this value. The resulting fit parameters for 
the remaining quantities are displayed in table HI Here, the quoted uncertainties include the 
statistical fit errors as well as a variation of Bkk in the range given by (j20p . From tabled] we 
infer that due to the reduction in the uncertainty of Bkk, correspondingly also the uncertainty 
in the slope parameters is much reduced, while the errors of the remaining parameters to a 
good approximation stay as before. Also the x^/n.d.f . is practically unchanged, still remaining 
below one for all values of Scut- 



■^Although the full Belle data set consists of 100 data points, we follow a suggestion of the experimentalists 
to only fit the data up to point 90 [48]. 
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Vs [GeV] 



Figure 1: Main fit result to the Belle data [24] for the differential decay distribution of the 
decay t~ — > Ut-Kstt". Our theoretical description corresponds to the fit of table H] with 
Scut = 4 GeV . The full fit including vector form factor F:^^{s) and scalar form factor 
Fq^'^{s) is displayed as the solid line. The separate vector and scalar contributions are shown 
as the dashed and dotted lines respectively. 



A graphical account of our central result is displayed in figured] The solid line corresponds 
to the fit of table m including vector form factor F^'^{s) as well as the scalar form factor 
Fq^{s) at Scut = 4 GeV . The separate contributions of F^^{s) and F^^{s) are shown as 
the dashed and dotted lines respectively. As is apparent, apart from the data points 5, 6 and 
7 in the low-energy region, our model provides a perfect description of the experimental data 
by the Belle collaboration. Also, from an inspection of the region below the K* resonance it is 
evident that a contribution from the scalar form factor Fq'^{s) is required, though, like in the 
analysis of ref. [36], the sensitivity to Fq '^(s) is not strong enough to allow for a determination 
of the corresponding model parameters. 



5 Conclusions 



Hadronic r decays provide a means to obtain information on low-energy QCD as well as 
hadron phenomenology. In this work, we have studied recent data on the decay channel 
T~ — > UrKsT^~ by the Belle collaboration [24]. The measured decay spectrum allows to 
test models for the vector and scalar Ktt form factors F^'"{s) and -Fg^'^(s), and to deduce 



the corresponding model parameters for the vector form factor. For F:^^{s), we have used 
a model which incorporates the constraints on the form factor from analyticity and elastic 
unitarity. Furthermore, we investigated n-subtracted dispersive integrals with a cutoff Scut, 
where n ranges from 1 to 4 and ^Scut was varied between 1.8 GeV and infinity. This allowed 
to test the sensitivity (or insensitivity) of the form factor model to contributions from higher 
energies which are not well known. For the scalar form factor Fq'^{s), we have employed 
the description of ref. [38], which is based on solving dispersion relations for a two-body 
coupled-channel problem, and was recently updated in [40]. 

Let us begin with summarising our final results for the parameters of the K* and K* 
vector resonances. As our central results, we quote the values of table [H at Scut = 4 GeV^. 
To the uncertainty given in table SI we add an error for the variation of our results when 
changing Scut- The resonance mass and width parameters are then found to be: 

niK* = 943.41 ± 0.59 MeV , 7x* = 66.72 ± 0.87 MeV , (21) 

ruj^,, = 1374 ± 45 MeV , 7^,/ = 240 ± 131 MeV , (22) 

while the mixing parameter for the second resonance reads 7 = — 0.039 it 0.020. As has been 



already stated above, the quantities of eqs. (pTj) and (j22|) are unphysical fit parameters. 

To obtain physical parameters for the K* and K* resonances, we have to compute the 
positions of the poles of the vector form factor in the complex s-plane. From the pole position 
Sp we can then read off the physical mass and width of the respective resonance according to 
the relation 

^ = Mr - ^ Tij . (23) 

Calculating the pole positions along the lines of the approach outlined in ref. [49] yields: 

Mii-. = 892.01 ± 0.92 MeV , F/^. = 46.20 ± 0.38 MeV , (24) 

Mj^,> = 1276 1^? MeV, F^,/ = 198 1^^ MeV. (25) 

The uncertainties are calculated by assuming a Gaussian error propagation while simultane- 
ously varying both m/j and 'Jr. The mass of the charged K* meson turns out rather close to 
the value advocated by the PDG, but more than 3 MeV lower than the Breit-Wigner reso- 
nance parameters obtained in refs. [24,36]. On the other hand, the K* width Tk* of eq. ()24p 
nicely agrees with the result of [24], but it is more than 4 MeV lower than the PDG average. 
To shed further light on the comparison with previous works let us calculate the pole 
position of the K*{892) for the best fit of ref. [36]. The model employed in this reference 
amounts, as far as the poles are concerned, to removing the term proportional to JieHxTris) 
from eq. (jlSp while keeping the energy dependent width. Denoting the respective fit param- 
eters with a tilde we have friK* = 895.28 ± 0.20 MeV and ^k* = 47.50 ± 0.41 MeV [36]. The 
corresponding pole position is given as the second line in table [5l being perfectly consistent 
with our results provided in eq. ()24p and the first line of table [5l The same exercise can be 
repeated for the original analysis performed by the Belle collaboration [24]. Here, we have 
employed the fit parameters corresponding to the second fit given in table 3, which are close 
to their final result for the K* parameters. In this case the corresponding pole position is 
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displayed as the last line in table [5j Again it turns out rather close to the previous results. 
To summarise, the pole position is found to be rather stable since different models yield com- 
patible values for the physical parameters Mk* and Tk* as defined in eq. (p3]) . Concerning 
the parameters of the second resonance, they are in general agreement with the findings of 
ref. [36], especially after the pole position is computed for the latter results. Due to the large 
uncertainties, however, we cannot make any more definite statement. 

Model Parameters Pole Positions 

{mK',7K*) [MeV] {Mk*,Tk*) [MeV] 

This work (943.41 ± 0.59, 66.72 ± 0.87) (892.0 ± 0.9,46.2 ± 0.4) 
Ref. [36] (895.28 ±0.20, 47.50 ±0.41) (892.1 ± 0.2,46.5 ± 0.4) 
Ref. [24] (895.47 ±0.20, 46.19 ±0.57) (892.5 ± 0.2,45.3 ± 0.5) 

Table 5: Comparison between model parameters and corresponding pole positions for the 
charged K*(892) meson. For definiteness, from ref. [24] we have employed the parameters of 
the second fit of their table 3 and consider only statistical uncertainties. 

Our fits to the r^ -^ VtKstt^ spectrum also allow to determine the slope parameters of the 
vector form factor F^'^{s). The advantage of using a three-subtracted dispersion relation is 
that the parameters A_|_ and A_|_ are directly determined from the data, making the extraction 
more model independent. The disadvantage being that therefore the uncertainties for \_^ 
turn out larger than for example in ref. [36] , where these parameters are a direct consequence 
of the form factor model. Higher slope parameters can of course also be calculated through 
dispersive integrals. For example in the case of A^ one has the relation: 

A'; = al + Za,a2 + M^ ^ Ids' %^ . (26) 

SK-n 

Together with the explicit fit results, this leads to 

A+ = (24.66 ± 0.77) • 10^^ A+ = (11.99 ± 0.20) • 10""^ , A+ = (8.73 ± 0.16) • 10"^ (27) 

where again the uncertainty due to the variation of Scut has been included in quadrature. 
Within the given errors, the value (j27p for A^ is in good agreement to the result of ref. [36], 
as well as the determination from an average of current experimental data for Ki^ decays [47]. 
On the other hand, both, the quadratic slope A_,_, and the cubic slope A^, are found somewhat 
lower than the corresponding results of ref. [36] . 

To conclude, differential decay spectra of hadronic r decays provide important information 
for testing form factor models and extracting the corresponding model parameters, thereby 
accessing QCD in the realm of low energies. It will be very interesting to see if our findings are 
corroborated by additional experimental data in the future. Furthermore, when comparing 
parameters of hadronic resonances, even when employing Breit-Wigner type parametrisations 
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with an energy-dependent width, pole positions in the complex s-plane should be provided in 
order to arrive at more model independent results. 
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